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Abstract 

With the discovery of the Higgs boson at the Large Hadron Collider the high energy physics 
community’s attention has now turned to understanding the properties of the Higgs boson, 
together with the hope of finding more scalars during run 2. In this work we discuss scenarios 
where using a combination of three decays, involving the 125 GeV Higgs boson, the Z boson 
and at least one more scalar, an indisputable signal of CP-violation arises. 

We use a complex two-Higgs doublet model as a reference model and present some benchmark 
points that have passed all current experimental and theoretical constraints, and that have cross 
sections large enough to be probed during run 2. 
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1 Introduction 


The discovery of the Higgs boson by the ATLAS [T] and CMS [ 2 ] collaborations at the Large 
Hadron Collider (LHC) has raised the interest of the high energy physics community in multi-scalar 
models. One of the most attractive features of some of these models is to provide extra sources of 
CP-violation which could help to explain the matter anti-matter asymmetry of the Universe. This 
was the reason that lead T.D. Lee to propose the two-Higgs double model (2HDM) [3] as a means 
to explain this asymmetry. Reviews of the 2HDM may be found, for example, in HE]- One of the 
CP-violating complex versions of the 2HDM, which we refer to as C2HDM, has been the subject 
of many recent studies [ 6 HI 2 ]. The C2HDM was first proposed in m and it is the simplest version 
of an explicit CP-violating 2HDM with a clear and easy limit leading to its CP-conserving version. 

As proposed in [T3], CP-violation in the scalar sector can be found in the interactions with gauge 
bosons in a very simple way. If CP were conserved, any decay hi hjZ would imply opposite CP 
parities for hj and hi. Moreover, assuming only lagrangian terms up to dimension four, any scalar hj 
decaying into ZZ would be CP even0. Thus, for example, the simultaneous presence of the decays 
hs —>• / 12 Z, /i 2 —>• hiZ, and h^ hiZ violates CP. We say that points in the C2HDM parameter 
space which lead to this situation belong to class Ci. Similarly (with the caveat in footnote 1), 
the simultaneous presence of the decays hi —)• hjZ, hi —?■ ZZ, and hj —>• ZZ, also violates CP. 
Within the 2HDM, there are three such possibilities, according to the {i,j) assignments, which we 
name classes C 2 , C 3 , and 6 * 4 . Notice that classes C'i-C '4 represent CP-violation, regardless of the 
origin of the neutral scalars. They may come from an N Higgs doublet model, or indeed from 
scalar fields in any number and from any representation of SU{2)l (singlets, doublets, triplets, 
combinations thereof, etc...) In Table [U we show the decays involved in each class. Furthermore, 


Classes 

Cl 

C2 

C3 

C4 

Cs 

Decays 

—>■ h2Z 
/i2 —y h\Z 
—y h\Z 

/i2 — y h\Z 
h\ —y ZZ 
/i2 —y ZZ 

—y h\Z 
hi —y ZZ 
/13 —y ZZ 

—>■ h2Z 
/i2 —y ZZ 
—y ZZ 

/13 —y ZZ 
/12 —y ZZ 
hi —y ZZ 


Table 1: Classes of combined measurements guaranteed to probe CP-violation in 2HDMs. 

in the specific context of a 2HDM, the properties of the fields ensure that, if CP were conserved, 
there would be two CP even neutral scalars and one CP odd neutral scalar, usually denoted by H, 
h, and A, respectively. Thus, in the 2HDM, the simultaneous presence of hi —?■ ZZ for i = 1,2,3 
signals CP-violation. We denote that possibility by class C 5 . We stress that class C 5 does not 
represent necessarily CP-violation in models other than the 2HDM. For example, even with three 
Higgs doublets one will surely have three neutral scalars and class C 5 would be consistent with 
CP-conservation. We will further discuss other classes that probe CP-violation that involve one 
scalar to two scalar decays that usually have the drawback of having smaller cross sections. 

It is interesting that there are only three basis-invariant quantities signalling CP-violation in the 
scalar sector of the 2HDM. They were introduced in [Slim, the connection with the observables 

^ There are CP conserving terms of dimension higher than four that can mediate the decay of a pseudoscalar 
into two vector bosons. Those could appear at loop level from a fundamental theory, but would lead to rates far 
smaller than the tree level rates considered in this article. A calculation performed in the framework of the 2HDM 
has shown m that the loop mediated decays of the type hi —^ ZZ are several orders of magnitude smaller than the 
tree-level ones. 
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explained in [T3], and revisited in [18]. Measurements of classes Ci-Cs are enough to probe all 
invariants. In the particular setting of the C2HDM, there is only one phase/source of CP-violation, 
all invariants are related, and the CP-violation in all classes (which one can take as the product of 
the three rates in each class) is proportional to that phase. 

One of the most interesting points of our proposal is that although the above described classes 
constitute an indisputable sign of CP-violation, they have all been searched for individually at run 
1. In fact, the searches hi ^ ZZ and hi —)• hjZ were already performed by both the ATLAS and 
CMS collaborations. Therefore, as long as we have enough signal events in three of the proposed 
channels for a given set of parameters, there are good chances of observing direct CP-violation at 
the next LHC run. 

This paper is organized as follows. In section [2l we briefly describe the complex 2HDM and 
the theoretical and phenomenological constraints imposed on the model with special emphasis on 
the most recent LHC data. In section [3l we propose a set of CP-violating benchmarks points for 
Type H and for the Flipped model. In the same section we discuss clear signs of CP-violation that 
involve the decay of one scalar to two scalars. Our conclusions are presented in Section 0] Finally, 
we present benchmark points for Type I and for the Lepton Specific model in appendix O 


2 The complex two-Higgs doublet model 


We use as a benchmark model an extension of the SM with an extra scalar doublet. This complex 
2HDM has a softly broken Z 2 symmetry (pi (pi,4>2 —4>2 and the scalar potential is written 

as [5] 


+ y + y + A 3 |(/>ip| 02 p -t A 4 {(t>\(t>2) 

+ + ( 1 ) 

and because the potential has to be hermitian, all couplings except m \2 and A 5 are real. In order to 
assure that the two phases cannot be removed simultaneously, we impose arg(A 5 ) 7 ^ 2 arg(mf 2 ) [13] . 
By taking m ^2 and A 5 real we recover the corresponding CP-conserving 2HDM. 

The model has three neutral particles with no definite CP, hi, /12 and h^, and two charged 
scalars The mass matrix of the neutral scalar states is obtained via the rotation matrix |19] 


R = 


C1C2 S1C2 S 2 ^ 

-(C1S2S3 -I-S1C3) C1C3 - S1S2S3 C2S3 
-C1S2C3 -I- S1S3 -(C1S3 S1S2C3) C2C3 


with. Si = sin Oi and Cj = cos ai {i = 1 , 2 ,3) and 


( 2 ) 


— tt/ 2 < ai < 7 r/ 2 , — 7 r /2 < 012 < '/r/2, —tt/2 < < vr/2. 


(3) 


The C2HDM has 9 independent parameters which we choose to be v, tan/3, ai, a 2 , 03 , 

mi, m 2 , and Re(m^ 2 )- With this choice the mass of heavier neutral scalar is a dependent parameter 
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given by 


m 


2 

3 


ml Ri3{Ri2 tan/3 - i?ii) + R23{R22 tan/3 - R 21 ) 

RasiRsi — R 32 tan/3) 


(4) 


and the parameter space will be restricted to values which obey m 3 > m 2 - 

We will analyse the usual four Yukawa versions of the C2HDM, in which the Z 2 symmetry is 
extended to the Yukawa Lagrangian [20] in order to avoid flavour changing neutral currents (FCNC). 
In all models the up-type quarks couple to (/>2 and the so-called Type I (Type II) is obtained by 
coupling down-type quarks and charged leptons to 4>2 while by coupling the down-type quarks 
to i/ii and the charged leptons to (j )2 we obtain the Flipped model and by coupling the down-type 
quarks to (/>2 and the charged leptons to (/)i we obtain the Lepton Specific model. 

We define the signal strength as 


hi _ (7 BR(/;.i —)■ /) 


where a is the Higgs boson production cross section and BR(/ij ^ /) is the branching ratio of the 
hi decay into the final state /; and BR®^(/i —>• /) are the corresponding quantities calculated 

in the SM. The cross sections were obtained from: HIGLU [21] - gluon fusion at NNLO, together 
with the expressions for the CP-violating model in |9|; SusHi [22| - bb ^ h at NNLO; [23| - Vh 
(associated production), tih and VV —)• h (vector boson fusion). The allowed parameter space of 
the C2HDM was recently reviewed in [TO] (see also p[^ll3lll9l[2TI[27] L The benchmark points that 
clearly signal CP-violation will be presented in the next section and are chosen from this set. The 
allowed points in parameter space are subject to the constraints we will briefly describe now. We 
note that we only focus here on scenarios where the lightest scalar hi is the 125 GeV Higgs. 


• We take the lightest neutral scalar, hi, to have a mass of 125 GeV in agreement with the 
latest results from ATLAS [28] and GMS [29] . 

• The accuracies in the measurements of the signal strengths in the processes pp —?■ hi —>• 
WW{ZZ), pp —/ii —77 and pp —/ii —)• r+T“ are about 20% at la [221130]. As shown 
in [0], imposing these run 1 constraints guarantees that the C2HDM automatically obeys all 
other run 1 constraints on the 125 GeV Higgs decays in this model. We will thus force 

and Pt-t to be within 20% of the expected SM value 

• The LHC results also allow us to put bounds on the heavier scalars /12 and / 13 . We impose 
the results on ^yy [31] in the range [145,1000] GeV and on [32] in the range [100,1000] 
GeV. We also use the results on hi —>■ ZZ —)• 4/ from [33] in the range [124,150] GeV and 
from [31] in the range [150,990] GeV, and on /i —>■ 77 from [341135] . Finally we also impose 
the constraints stemming from the results based on the searches hi —?■ Zhi —>• Zbb{T~^T~) [36] 
and hi —>• Zhi —)■ llbb m- 

• We consider the constraints on the charged Higgs Yukawa vertices that depend only on the 
charged Higgs mass and on tan/3. There is a new bound on 6 —?• 57 , in Type H/F [33] of 
'mu± > 480 GeV at 95% C.L.. Putting together all the constraints from B-physics [321140] 
and also from the Rb = T{Z —)• bb)/T{Z —)• hadrons) [41] measurement, we can state that 
roughly tan /3 > 1 for all models. LEP searches on e'’'e“ —)• H^H~ [42] and the LHC searches 
on pp —)• it{-^ H^b [431144] 1 lead us to roughly consider mjj± > 100 GeV in Type I/LS. 
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• We consider the following theoretical constraints: the potential has to be bounded from 
below [l5], perturbative unitarity is required [IHl - HH] and all allowed points comply with the 
oblique radiative parameters [iSHST] . 

• The scenarios we will present in the next section are a clear signal of CP-violation in models 
with an extended scalar sector. Models with a CP-violating scalar sector are constrained 
by bounds from electric dipole moments (EDMs) measurements. Although the search for 
the proposed final states should be performed from a model independent perspective, we 
will nevertheless estimate the most important constraints on the CP-violating phases in the 
context of the C2HDM [71 I52H56] . 

The most stringent bound [7] comes from the ACME 15Z1 results on the ThO molecule EDM. 
In order to have points with EDMs of an order of magnitude that conforms to the ACME 
result, we have computed the Barr-Zee diagrams with fermions in the loop. As we will see, the 
ACME bound can only be evaded by either going to the limit of the CP-conserving model or 
in scenarios where cancellations [55ll56j among the neutral scalars occur. These cancellations 
are due to orthogonality of the R matrix in the case of almost degenerate scalars [S]. We 
should hnally point out that ref. [53] argues that the extraction of the electron EDM from 
the data is filled with uncertainties and an order of magnitude larger EDM than that claimed 
by ACME should be allowed for. 

3 CP-violating benchmark points 

In this section, we present some benchmark points that allow us to definitely probe CP-violation 
during LHC’s run 2. In table [2l we present four benchmark points, where the hrst three are for 
Type II and PA is for the Flipped model (Type I and Lepton Specific are discussed in appendix A). 
For each point we give the values of the parameters of the model, the values of the pseudoscalar 
component of the Yukawa coupling of the lightest Higgs and the values of the cross sections for 
the different processes. The cross sections are calculated assuming that all scalars in the final 
state are detected in the decay to bb and all Z bosons are detected in the leptonic decays, providing 
therefore a very conservative estimate for the number of signal events available. Regarding the cross 
sections, we sum over all possible production process with one scalar in the final state. Therefore, 
the numbers presented in the table correspond either to 

a{pp hi + X) BR{hi —>• hjZ) BR{hj —>• bb) BR{Z —)■ ll) , (6) 

or 

a{pp^ hi + X)BR{hi^ ZZ)BR^{Z ^ll) (7) 

and I = e, p. 

The general criteria for the choice of our benchmark points is the following: the points have 
passed all the constraints described in the previous section; the number of events for a luminosity 
of 100/6“^ should be at least above 50, and the smallest number in table [2] for this luminosity is 
61 events. Note that this number already takes into account the decay of the scalar into bb and the 
decay of all Z bosons into leptons (a reduction of 0.06732 for each Z). Therefore, we expect a much 
larger number of events when all other combinations of final states are taken into account by the 
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PI 

P2 

P3 

P4 

ai 

1.12569 

1.04842 

-1.33589 

1.41610 

02 

0.49091 

-0.00825 

-0.00129 

0.24037 

03 

-1.56775 

0.00674 

0.63749 

-0.81993 

/3 

0.92913 

1.00182 

1.27669 

1.29413 

tan /3 

1.33845 

1.56366 

3.30155 

3.52171 

mi (GeV) 

125.00 

125.00 

125.00 

125.00 

m 2 (GeV) 

127.32 

273.15 

282.53 

231.74 

m3 (GeV) 

252.63 

421.64 

287.80 

360.59 

mH± (GeV) 

481.25 

452.50 

604.89 

527.67 

Re(m?2) (GeV)2 

-0.5625E-k02 

0.1183E-k05 

0.1590E-k05 

0.2156E-k05 

boi 

-0.63099 

0.01291 

0.00426 

-0.83837 

bhi 

-0.63099 

0.01291 

0.00426 

0.06760 

Ci[l] cr3xBR(h3 —>■ /12Z — bbll) 

114.528 [fb] 

61.529 [fb] 

0.000 [fb] 

27.484 [fb] 

Cl [2] cr2xBR(/i2 -^hiZ ^ bill) 

0.000 [fb] 

0.615 [fb] 

7.401 [fb] 

18.462 [fb] 

Cl[3] (T3xBR(/i3 — >• hiZ — bbll) 

26.656 [fb] 

1.100 [fb] 

24.519 [fb] 

1.787 [fb] 

C2[l] (T2xBR(h2 — >■ hiZ — ^ bbll) 

0.000 [fb] 

0.615 [fb] 

7.401 [fb] 

18.462 [fb] 

C2 [2] CTi X BR(hi ^ ZZ ^ nil) 

5.495 [fb] 

5.792 [fb] 

5.592 [fb] 

4.802 [fb] 

C2[3] cr2xBR(h2 ^ ZZ ^ nil) 

1.386 [fb] 

2.598 [fb] 

1.802 [fb] 

1.220 [fb] 

C3[l] cr 3 xBR(/i 3 — >■ hiZ — bbll) 

26.656 [fb] 

1.100 [fb] 

24.519 [fb] 

1.787 [fb] 

C3 [2] CTi X BR(hi ^ ZZ ^ nil) 

5.495 [fb] 

5.792 [fb] 

5.592 [fb] 

4.802 [fb] 

C3[3] cr3xBR(/i3 ^ ZZ ^ nil) 

1.011 [fb] 

0.003 [fb] 

1.733 [fb] 

1.058 [fb] 

C4[l] cr3xBR(h3 —>■ h 2 Z — bbll) 

114.528 [fb] 

61.529 [fb] 

0.000 [fb] 

27.484 [fb] 

C 4 [ 2 ] (72 xBR(h2 ^ ZZ ^ nil) 

1.386 [fb] 

2.598 [fb] 

1.802 [fb] 

1.220 [fb] 

C4[3] (T3xBR(/i3 ^ ZZ ^ nil) 

1.011 [fb] 

0.003 [fb] 

1.733 [fb] 

1.058 [fb] 

C 5 [l] (T3xBR(/i3 ^ ZZ ^ nil) 

1.011 [fb] 

0.003 [fb] 

1.733 [fb] 

1.058 [fb] 

C 5 [ 2 ] (72xBR(h2 ^ ZZ ^ nil) 

1.386 [fb] 

2.598 [fb] 

1.802 [fb] 

1.220 [fb] 

C 5 I 3 ] aixBR{hi ^ ZZ ^ nil) 

5.495 [fb] 

5.792 [fb] 

5.592 [fb] 

4.802 [fb] 


Table 2: Benchmark points for Type II: PI, P2 and P3, and for the Flipped model: P4, for LHC at ^/s = 13 
TeV. All Z bosons decay leptonically which corresponds to a factor of 0.06732 for each Z decay. 


experiments (as it is obviously the case for the ZZ hnal states, where we can have combinations 
of leptons and jets final states). In table [3] we show the rates obtained for the benchmark points 
which are then compared to the available experimental data from the LHC at y/s = 8 TeV. 

The criteria for the choice of each particular point is severely constrained by the ACME results. 
In fact, all the points have similar features in that they either have two neutral scalar masses 
almost degenerate or values of the angles very close to zero (therefore approaching the limit of the 
CP-conserving 2HDM). Points PI and P3 have degenerate masses while point P2 has very small 
02 and 03 values. That is why for point P2, the decay /12 —>■ hiZ is suppressed. In the limit 
02 = 03 = 0, /13 is the pseudo-scalar and hi and h 2 are scalars and h 2 ^ hi Z is forbidden. For the 
same reason, h^ ^ ZZ is forbidden. Note however that although 02 and 03 are very small we still 
have a large number of signal events for 100/5“^ in /12 —)• hiZ. As 02,3 move away from zero (the 
CP-conserving limit) certain CP-violating observables grow extremely fast. Thus, we can be very 
close to this limit and still have large CP-violating signals. 

The points were also chosen so that they would probe more than one class simultaneously. PI 
probes classes C 3 , 6*4 and C 5 ; P2 probes Ci and C 2 ; P3 probes C 2 , C 3 and C 5 while the point for 
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PI 

P2 

P3 

P4 

^J,ww{hl) = ^izz{hl) 

1.09016 

1.14962 

1.11696 

0.95402 

l^TT (^1 ) 

1.16717 

0.98826 

0.96621 

1.02628 

M 77 (^ 1 ) 

0.92139 

1.02589 

0.87922 

0.85345 

^J'bb(VH){hl) 

0.71662 

0.93593 

0.65922 

0.94294 

mz(h2)liiTz 

0.225/NA 

0.225/1.264 

0.151/0.185 

0.151/0.190 

0.117/0.170 

0.117/0.176 

0.058/0.121 

0.058/0.130 

[lrr(h2)l\lTr 

1.59/ 3.98 

180.00/ 472.37 

7.98/ 490.42 

0.90/ 363.88 

aBR^^{h2)/aBR^^:P [fb] 

15.265/ 29.705 

0.318/2.678 

0.011/2.727 

0.018/5.998 

Ai^^(/i 2 )/Ai“P (to < 150GeV) 

aBRzh^zbb{h2)/ [pb] 

0.258/0.259 
0 . 000 / 0.000 

0 . 000 / 0.000 

0.003/0.308 

0 . 000 / 0.000 

0.042/0.250 

0 . 000 / 0.000 

0.108/0.403 

aBRzh^Zrrih2)/aBR^^^^^^^ [pb] 

0 . 000 / 0.000 

0.000/0.105 

0.005/0.089 

0.012/0.085 

aBRzh^iibb{h2)/o'BR^^^Hf^^ [fb] 

0 . 000 / 0.000 

0.222/ 15.242 

2.855/ 12.167 

7.259/ 14.082 

fJ-ww {h^) / 

mz{h^)li^Tz 

0.053/0.074 

0.053/0.068 

0.000/0.083 

0.000/0.086 

0.111/0.125 

0.111/0.147 

0.072/0.099 

0.072/0.095 

Mrr(/l3)/Mrr'' 

3.12/ 427.59 

8.70/ 1241.83 

13.52/ 500.43 

0.04/ 663.64 

aBR^^{h3)laBR^^^P [fb] 

0.022/6.511 

0.028/2.002 

0.010/2.672 

0.004/2.823 

aBRzh^zbb{h3)/aBR‘'zh^zbb [pb] 

aBRzh^ZTT{h3)/aBR^^P^^^^ [pb] 

0.147/0.310 

0.017/0.102 

0.005/0.081 

0.001/0.035 

0.135/0.228 

0.016/0.081 

0.009/0.156 

0.001/0.038 

aBRzh-yiibb{h3)/o'BR^I^^ii^f^ [fb] 

9.926/ 23.839 

0.337/2.731 

9.076/ 15.230 

0.605/7.358 


Table 3: Constraints from the LHC at y/s = 8 TeV for the benchmark points PI, P2 and P3 (Type II) and 
P4 (Flipped). NA stands for not available. 


the Flipped model probes all classes. Furthermore, points PI and P4 were also chosen to show 
that large pseudoscalar components are not only still allowed, as previously discussed in m, but 
they can also easily be probed at the next LHC run. 

Finally, in table[l]we present the production cross sections for hi, /12 and h^. In the same table 
we show the (j{hi) x Br{hi —)■ X) where X stands for the main final sates being searched by ATLAS 
and CMS at the next LHC run. These numbers allow the experimental groups to understand if a 
given scalar is found in direct production whether it comes from a CP-violating process or not. In 
the same table we also present the values of the scalar production cross sections that lead to decays 
of the type hi —)• hjhj and hi —)• hjhk and that are clearly too small to be detected at the LHC for 
the sets of benchmark chosen, except for a few cases for points PI and P4. 

3.1 CP-violating scenarios involving scalar to two scalars decays 

There are other classes of decays that constitute a sign of CP-violation in the 2HDM. Some of them 
involve the decay h^ —)• h 2 hi which is not present in the CP-conserving version of the 2HDM. In 
fact the decay h^ —)• h 2 hi is only possible if either all hi are CP-even, two of the hi are CP-odd 
and one is CP-even or if CP is not conserved. Since the decay hi —)• ZZ was already observed 
we know hi has a CP-even component. Therefore, we can discuss the combinations of decays that 
together with h^ —)• h 2 hi will be a clear sign of CP-violation in the 2HDM or that will point to 
other extensions of the SM that can be either CP-conserving or CP-violating. 

In table [5] we present new classes of decays that again constitute model independent signs of 
CP-violation. Class Cq is composed by the three decays hi —)• ZZ, h^ / 12/11 and /i 3 —>■ h 2 Z. 
There are other sets of simultaneous measurements involving /13 ^ h 2 hi that are consistent with 
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PI 

P2 

P3 

PA 

(j(/ii) 13TeV 

61.600 

[pb] 

53.217 

[pb] 

54.825 

[pb] 

51.275 

[pb] 

cr(/ii)BR(/ii ^ W*W*) 

11.819 

[Pb] 

12.459 

[pb] 

12.028 

[pb] 

10.328 

[pb] 

a{hi)BR{hi Z*Z*) 

1.212 

[pb] 

1.278 

[pb] 

1.234 

[pb] 

1.060 

[pb] 

a{hi)BR{hi —>• bb) 

34.383 

[pb] 

29.087 

[pb] 

28.256 

[pb] 

30.313 

[pb] 

a{hi)BR{hi —>■ tt) 

3.969 

[pb] 

3.360 

[pb] 

3.264 

[pb] 

3.485 

[pb] 

(t(/ii)BR(/ii —>■ 77) 

129.973 [fb] 

144.664 [fb] 

123.188 [fb] 

120.222 [fb] 

(72 = <7(112) 13TeV 

56.583 

[pb] 

4.262 

[pb] 

1.602 

[pb] 

3.354 

[pb] 

Cr2XBR(/l2 ^ WW) 

2.814 

[pb] 

1.323 

[pb] 

0.910 

[pb] 

0.656 

[pb] 

(72XBR(/l2 —>■ ZZ) 

0.306 

[pb] 

0.573 

[pb] 

0.398 

[pb] 

0.269 

[pb] 

(T2xBR(Ii2 —>■ bb) 

42.534 

[pb] 

1.894 

[pb] 

0.067 

[pb] 

1.944 

[pb] 

(72 xBR(ft,2 —>■ tt) 

4.911 

[pb] 

0.224 

[pb] 

0.008 

[pb] 

0.002 

[pb] 

(T2XBR(Il2 —>■ 77) 

35.041 

[ib] 

0.879 

[fb] 

0.027 

[ih] 

0.046 

[ih] 

(T2xBR(Ii2 —>■ hiZ) 

0.000 

[pb] 

0.017 

[pb] 

0.213 

[pb] 

0.464 

[pb] 

fT2xBR(/i2 —>■ hiZ —>■ 66Z) 

0.000 

[pb] 

0.009 

[pb] 

0.110 

[pb] 

0.274 

[pb] 

(T2xBR(/i2 —>■ hiZ —>■ ttZ) 

0.000 

[fb] 

1.055 

[fb] 

12.697 

[fb] 

31.530 

[fb] 

(T2xBR(/i2 —>■ hihi) 

0.000 

[fb] 

0.007 

[fb] 

5.016 

[fb] 

0.000 

[fb] 

(T2xBR(/i2 —>■ ^1^1 —>■ bb bb) 

0.000 

[fb] 

0.002 

[fb] 

1.332 

[fb] 

0.000 

[fb] 

(j2xBR(/i2 —>■ /ll/ll —>■ bb tt) 

0.000 

[fb] 

0.000 

[fb] 

0.308 

[ih] 

0.000 

[ih] 

(72xBR(/i2 —^ hihi —>■ TT tt) 

0.000 

[fb] 

0.000 

[fb] 

0.018 

[fb] 

0.000 

[fb] 

(J3 = ct(/i3) 13TeV 

4.043 

[pb] 

8.480 

[pb] 

2.086 

[pb] 

1.819 

[pb] 

Cr3XBR(/l3 ^ ITIT) 

0.526 

[pb] 

0.001 

[pb] 

0.871 

[pb] 

0.509 

[pb] 

(j3xBR(/i3 —>• ZZ) 

0.223 

[pb] 

0.001 

[pb] 

0.382 

[pb] 

0.233 

[pb] 

(T3xBR(Ii3 —>• bb) 

0.047 

[pb] 

0.016 

[pb] 

0.109 

[pb] 

0.058 

[pb] 

(T3 xBR(ft,3 —>■ tt) 

5.558 

[fb] 

1.913 

[fb] 

12.856 

[fb] 

0.020 

[fb] 

(T3xBR(/i3 77) 

0.059 

[fb] 

0.093 

[fb] 

0.028 

[fb] 

0.013 

[fb] 

(T3xBR(Ii3 —>• hiZ) 

0.709 

[pb] 

0.030 

[pb] 

0.707 

[pb] 

0.045 

[pb] 

(T3xBR(/i3 —>• hiZ —>• bbZ) 

0.396 

[pb] 

0.016 

[pb] 

0.364 

[pb] 

0.027 

[pb] 

(T3xBR(/i3 —>■ hiZ —>■ ttZ) 

45.708 

[fb] 

1.887 

[fb] 

42.067 

[fb] 

3.051 

[fb] 

(T3xBR(/i3 —>• h 2 Z) 

2.263 

[pb] 

2.057 

[pb] 

0.000 

[pb] 

0.705 

[pb] 

(T3xBR(/i3 —>• h 2 Z —>• bbZ) 

1.701 

[pb] 

0.914 

[pb] 

0.000 

[pb] 

0.408 

[pb] 

(T3xBR(/i3 —>■ h 2 Z —>■ ttZ) 

196.416 [fb] 

107.996 [fb] 

0.000 

[fb] 

0.500 

[fb] 

(T3xBR(/i3 —>• hihi) 

0.090 

[fb] 

0.230 

[fb] 

2.071 

[fb] 

19.918 

[fb] 

(T3xBR(/i3 hihi —>■ bb bb) 

0.028 

[fb] 

0.069 

[fb] 

0.550 

[fb] 

6.961 

[fb] 

(j3xBR(/i3 —>• hihi —>■ bb tt) 

0.007 

[fb] 

0.016 

[fb] 

0.127 

[fb] 

1.601 

[fb] 

(73XBR(/l3 —>■ hihi —>■ TT tt) 

0.000 

[fb] 

0.001 

[fb] 

0.007 

[fb] 

0.092 

[fb] 

(T3XBR(/l3 —>• /l2^l) 

263.916 [fb] 

0.038 

[fb] 

0.000 

[fb] 

11.157 

[fb] 

(T3xBR(/i3 —>• /12/11 —>■ bb bb) 

110.732 [fb] 

0.009 

[fb] 

0.000 

[fb] 

3.822 

[fb] 

(j3xBR(/i3 —>• /12/11 —>■ bb tt) 

25.567 

[ih] 

0.002 

[fb] 

0.000 

[ih] 

0.444 

[ih] 

(73xBR(/i3 —>■ h 2 hi —>■ TT tt) 

1.476 

[fb] 

0.000 

[fb] 

0.000 

[fb] 

0.001 

[fb] 


Table 4: Predictions for a x BR for the LHC at ^/s = 13 TeV for the benchmark points PI, P2 and P3 
(Type II) and P4 (Flipped). 


CP conservation in models with more than two Higgs doublets, and which allow the determination 
of the possible CP assignments. These are; 

• [hi —)■ ZZ; /13 —)■ h 2 hi] /13 —)■ hiZ^ leading to the CP assignments (+, —, —); 
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• \hi ZZ\ /13 —)• h 2 hi] /i 2 —>• hiZ, leading to the CP assignments (+, —, —); 

• \hi —)• ZZ-, /13 —>■ h 2 hi] h 2 —)• ZZ, leading to the CP assignments (+, +, +). 

There are still other combinations involving scalar to scalar decays that are model independent 
signs of CP-violation. Such is the case of class C^ in table [5] composed by the decays /i 2,3 —)• hi hi, 
h 2,3 —^ hiZ and hi —)• ZZ. Finally, other combinations like hs —>■ hihi(h 2 h 2 ), h 2 —)• hi hi and 
hi —)• ZZ are not possible in a CP-conserving 2HDM but are possible in the C2HDM and can also 
serve to determine the CP-quantum numbers of other extensions of the scalar sector. A detailed 
study of these classes will be performed in a forthcoming publication [58]. 


Classes 

Ce 

C7 

Decays 

/13 —>■ /12/11 
/13 —>■ /l2^ 
hi —y ZZ 

^2,3 hihi 

^2,3 hiZ 

hi —y ZZ 


Table 5 : Classes of combined measurements guaranteed to probe CP-violation. 


4 Conclusions 

We have proposed hve classes of processes that constitute conclusive evidence of CP-violation in 
scalar decays. While the C 5 class is particular to the C2HDM, all other classes are valid in any 
scalar extension of the SM. One of the most attractive features of our proposal is to rely on searches 
that are already planned for the LHC run 2, namely h* —)• ZZ and hi —)• hjZ. Furthermore, it 
does not depend on complex distributions nor asymmetries of any kind, but only on total rates of 
specific processes. It is a direct and straightforward way to search for CP-violation at the LHC in 
scalar decays. As far as we know this is the only method of probing CP-violation based on rates 
only. 

We have shown that even taking into account all constraints and in particular the one from 
ACME that heavily restricts the amount of CP-violation in the model, it is still easy to find points 
to probe each of the proposed classes. In many cases a point can be used to probe several classes 
simultaneously. We have chosen a set of benchmark points according to different criteria, always 
keeping in mind that the decays should be within the reach of the LHC’s run 2. We should point out 
however that even if these points are excluded the parameter space is large enough to provide many 
more points and the model is far from excluded (nor is CP-violation in scalar decays excluded). 
The future bounds on EDMs [521159] can have a strong impact on the allowed parameter space, and 
one has to consider the interplay between the EDM bounds and the data from run 2 to propose 
scenarios for future experiments. However, the EDM constraints get looser if one goes beyond 
the setting discussed here, allowing for Ag / 0 and/or A 7 7 ^ 0 [11]. In that case classes Ci to C 5 
still probe CP-violation, and thus the methods proposed here should be pursued experimentally 
regardless of the fate of the C2HDM. In particular, classes Ci to C 4 probe CP-violation in all 
models. 

We also propose two new classes of decays, Cg and C 7 that involve the already observed decay 
hi — ZZ, one decay of the type hi —)• hjhj(j^-^ with j ^ k and one decay of type hj —)• hkZ. 

As important guidelines for experiments, we propose six benchmark points covering all C2HDM 
types: type H (P1-P3), and Flipped (P4) in tables [21 [3] and m type I (P5) and Lepton Specific 






(P6) in tables [H [7] and [8] of appendix A. We provide all event rates for all scalar processes and for 
each benchmark points. This allows not only to search for the CP-violating classes of decays but 
also to confirm or disprove the points via direct search of each scalar. If a particular point is found 
the other decays could clarify if we are in presence of the C2HDM or of some other CP-violating 
extension of the SM. 


A Benchmark points for Type I and for the Lepton Specific models 



P5 

P6 

ai 

1.30680 

1.08742 

Oi2 

0.10867 

0.00960 

03 

-0.20624 

-0.41962 

P 

1.15333 

1.03051 

tan P 

2.25459 

1.66717 

mi (GeV) 

125.00 

125.00 

m 2 (GeV) 

235.45 

262.98 

m3 (GeV) 

359.20 

264.60 

mH± (GeV) 

522.87 

471.76 

Re(m22) (GeV)2 

0.9504E-^02 

-0.3006E-H05 

boi 

0.04810 

0.00576 

bLi 

0.04810 

-0.01600 

Ci[l] (T3xBR(h3 —>■ / 12 Z —>• bbll) 

1.251 [fb] 

0.000 [fb] 

Ci[2] (T2xBR(h2 —>■ hiZ —>• bbll) 

5.644 [fb] 

3.030 [fb] 

Ci[3] (T3xBR(/i3 —>• hiZ —>• bbll) 

15.477 [fb] 

27.984 [fb] 

C2[l] (T2xBR(h2 — >■ hiZ — >• bbll) 

5.644 [fb] 

3.030 [fb] 

C 2 [2] aixBR{hi^ ZZ ^ llll) 

4.954 [fb] 

5.146 [fb] 

C2[3] cr2xBR(/i2 ZZ ^ nil) 

1.934 [fb] 

1.053 [fb] 

C3[l] (T3xBR(h3 — >• hiZ — >• bbll) 

15.477 [fb] 

27.984 [fb] 

C3[2] aixBR{hi ^ ZZ ^ llll) 

4.954 [fb] 

5.146 [fb] 

C3[3] cr3xBR(/i3 ZZ ^ nil) 

1.326 [fb] 

1.840 [fb] 

C4[l] (T3xBR(h3 — >• h2Z — >• bbll) 

1.251 [fb] 

0.000 [fb] 

Ca[ 2 ] (T 2 xBR(/i 2 zz ^ nil) 

1.934 [fb] 

1.053 [fb] 

C4[3] Cr3XBR(/l3 ZZ ^ nil) 

1.326 [fb] 

1.840 [fb] 

Cs)!] Cr3XBR(/l3 ZZ ^ nil) 

1.326 [fb] 

1.840 [fb] 

C 5 [ 2 ] (j 2 xBR{h 2 ZZ ^ nil) 

1.934 [fb] 

1.053 [fb] 

C 5 [3] aixBR{hi^ ZZ ^ nil) 

4.954 [fb] 

5.146 [fb] 


Table 6: Benchmark points for Type I: P5 and for the Lepton Specific model: P6, for LHC at ^/s = 13 TeV. 
All Z decay leptonically corresponding to a factor of 0.06732. 

In this appendix we present two further benchmark points, one for Type I and the other for 
the Lepton Specific (LS) model. In table [6] we present the values of the parameters and the cross 
sections for benchmark point P5 in Type I and P6 for the LS model. In Type I it was possible to 
find a point that not only complies with all the constraints but that probes all CP-violating classes 
at the same time. For the LS model the classes probed are C 2 , C 3 and C 5 . 

As we did for the remaining benchmark points, we present in table [7] the effect of the LHC 
constraints on the processes involving scalars. In table [8] we present the production cross sections 
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P5 

P 6 

gww{hi) = pzzihi) 

0.98240 

1.02070 

fJi-rr (^1) 

1.12419 

0.83628 

/^77 (^1) 

0.84875 

0.86872 

gbb{VH){hl) 

0.99480 

1.02881 

gww{h2)/ 

gzz{hF)lpTz 

0.091/0.115 

0.091/0.111 

0.058/0.108 

0.058/0.112 

pirr{h2)lgV7 

0.56/ 377.80 

72.97/ 451.42 

aBR.y^{h2)/aBR^^P [fb] 

0.046/3.975 

0.125/6.838 

g.y-,{h 2 )/g^^P (m < 150GeV) 
aBRzh^Zbb{h2)/crBRfP^Zbb [pb] 

0 . 000 / 0.000 

0.032/0.337 

0 . 000 / 0.000 

0.016/0.349 

(TBRzh^ZTT{h2)lcFBR^^l^^^^ [pb] 

0.004/0.080 

0.001/0.114 

aBRzh^iibb{h2)/ [fb] 

2.127/ 13.013 

1.100/ 27.341 

gww {h'i) / 

gzz{h^)lgTz 

0.087/0.097 

0.087/0.094 

0.102/0.113 

0.102/0.123 

grr{h:i)lgTP 

0.89/ 656.23 

281.79/ 454.89 

aBR.y^{h3)/aBR^^P [fb] 
aBRzh^zbb{h3)/aBR^Zh^zbb [pb] 
CrRi?z/i^Zrr(h3)/crRi?|’)f_^2rr [pb] 

aBRzh^iibb{h3)/aBR2;^_^iiiji, [fb] 

0.046/2.758 

0.075/0.155 

0.009/0.038 

5.077/7.483 

0.875/6.334 
0.151/0.348 
0.013/0.114 
10.163/ 24.919 


Table 7: Constraints from the LHC at y/s = 8 TeV for the benchmark points P5 (Type I) and P6 (Lepton 
Specific). 


for hi, /i 2 and and also the a{hi) x Br{hi —)• X) where again X stands for the most relevant final 
states searched by ATLAS and CMS at the next LHC run. We also show the values of the scalar 
production cross sections that lead to decays of the type hi hjhj and hi —>■ hjh^. Interestingly, 
for the benchmark points of Type I and Lepton Specific, there are many scalar to scalar decays 
that could be probed at the next LHC run. 
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